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Abstract. - We discuss a recent experiment in which the spectrum of the vortex line density 
fluctuations has been measured in superfluid turbulence. The observed frequency dependence 
of the spectrum, f~^^^, disagrees with classical vorticity spectra if, following the literature, the 
vortex line density is interpreted as a measure of the vorticity or enstrophy. We argue that the 
disagrement is solved if the vortex line density field is decomposed into a polarised field (which 
carries most of the energy) and an isotropic field (which is responsible for the spectrum) . 



^ ', Motivation and aim. — Recent experiments have 
Q explored similarities and differences between turbulence 
Qjin classical, ordinary fluids and turbulence in He II (su- 
' 'perfluid turbulence). Superfluid turbulence consists of a 
I tangle of quantised vortex filaments; it is usually charac- 
^terised (in both experiments and numerical simulations) 
00 by the vortex line density L (defined as the vortex length 
per unit volume) . Superfluid turbulence can be generated 
^Din many ways: heat currents [1-3], vibrating wires [4], 
• oscillating grids [5] or spheres [6], towed grids [7], bel- 
^^lows [8,9], rotating propellers [10,11] and ultrasound [12]. 

Superfluid turbulence is also tackled in the slightly differ- 
^^ent context of superfluid •^He-B [13,14]. 

^ \ The current understanding of superfluid turbulence at 
k>l relatively high temperature is the following. According 
rStiO experimental [10,15], theoretical [16] and numerical [17] 
results, at sufHciently large scales in the inertial range, the 
liormal fluid and the superfluid components of He II are 
strongly coupled, the superfluid and normal fluid velocities 
are matched, ~ w„, and their energy spectra obey the 
classical Kolmogorov law k~^^^ (where the wavenumber k 
and the frequency / are related by = f /V where V is 
the mean flow). 

In a recent experiment, Roche et al. [11] measured the 
spectrum of the fluctuations of the vortex line density L in 
turbulent superfluid helium "^He at T = 1.6 K and found 
a clear /^^^'^ dependence. Our aim is to reconcile this 
observation with the current understanding of superfluid 
turbulence and the interpretation (which is quite common 
in the literature) of the vortex line density L as a measure 
of the superfluid vorticity, uig — kL, where k « 10~^m^/s 
is the quantum of circulation. Under this interpretation, 
the /"^/"^ spectrum observed by Roche et al. [11] seems 



to contradict the scaling of vorticity observed in classical 
turbulence, which is a flat or slowly decreasing frequency 
spectrum (see for example [18,19] and references within). 

To be more precise, what was actually measured in 
Ref. [11] is the vortex line density corrected by a sine 
squared prefactor to account for the orientation of each 
vortex line (as explained in Ref. [3] for example) : the com- 
ponent of a vortex line parallel to the direction of sound 
propagation does not contribute to the second sound at- 
tenuation at flrst order. The natural interpretation of the 
measured signal is therefore the magnitude of the vorticity 
corrected by a prefactor calculated with the orientation of 
the vorticity vector. Using the DNS dataset"'^ of Gotoh et 
al. [20], we have checked that the correction introduced by 
this prefactor has only a small contribution to the slope of 
the spectrum of the magnitude of the vorticity in classical 
turbulence. This still leaves us with a major discrepancy 
between the classical vorticity spectrum and the steeper 
j-5/3 decrease which was observed. 

In bringing together the current understanding of super- 
fluid turbulence with the observed spectrum of L, we also 
need to make sure that the interpretation of all measure- 
ments performed in Ref. [11] using a pressure sensor and 
a second sound detector are consistent with each other. 
In fact, at first sight there seems to be an inconsistency 
between the mean vortex line density and the energy of 
the flow estimated from the measured velocity 

For a mean velocity V « Im/s at T = 1.6 K, Roche et 
al. [11] report a mean vortex line density L corresponding 
to an average intervortex spacing S 



^Data were downloaded from International Computational Fluid 
Dynamics database, hosted by the Cineca supercomputing center, 
Bologna, Italy, (http://cfd.cineca.it/) 
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(5= l/\/Z«4x IQ-^m. (1) 
from which we estimate backwards: 

1= 1/^2 «6 X 10^°m-2. (2) 

The kinetic energy per unit volume of the same flow is 

K = \pnvl + \psvl. (3) 

If the normal fluid and the supcrfluid arc indeed coupled, 
Vn ~ Vs, and if we assume that this velocity is approxi- 
mately Vrms = 0.3m s~^ (corresponding to the mean flow 
V = Im s~^ and a turbulence intensity of 30%), we have 

K « \{Pn + Ps)V,^ms = ^Vrms = S-SJ m'^. (4) 

where p = ps + Pn = 145 kg m"''. Since pn < 0.2ps, this 
energy K is approximately equal to the superfluid kinetic 
energy. Let us assume that the kinetic energy per unit 
volume, is approximately equal to the kinetic energy 
per unit length, £, times the length per unit volume L : 

K w £L. (5) 

The kinetic energy per unit length is obtained in cylindri- 
cal coordinates (r, 0, z) by integrating the square of the ve- 
locity field K/(27rr) around a straight vortex line (set along 
the z direction) from the radial distance r = a k, 10~^'^m 
(the vortex core radius) to some upper cutofi: b: 

If we take b = S, using ps = 122 kg m~^ at T = 1.6K, we 
have 

f » 1.0 X 10-^^J m-\ (7) 

thus 

L » K/£ = 6 X (8) 

which is much bigger than the value of L from second 
sound measurements, L = 6 x lO^^m"^. This second ap- 
parent inconsistency must be solved too. 

The following model which we propose to solve these 
inconsistencies may not be the final answer. Nevertheless, 
we think that the exercise of putting together a coher- 
ent scenario with the experimental information which is 
available at this stage is a valuable exercise which should 
stimulate further work and bring us closer to the correct 
solution of the puzzle. 

Model. — The model which wc propose to solve the 
puzzle described in the previous section has two key 
features: the decomposition of the vortex line distri- 
bution into a "polarised" field and a "unpolariscd" (or 
"isotropic" ) field, and the assumption that the unpolarised 
field has some statistical properties of a passive vector 
field. We stress that our interpretation is preliminary. 



Decomposition of the vortex line density. Since the 
vortex core radius is many orders of magnitude smaller 
than 5 or any other length scale of interest in the flow, 
we follow Schwarz [21] and describe vortex lines as space 
curves s(^, t) where ^ is arclcngtli and t is time. The quan- 
tity s' = ds/d^ is the unit vector at the point x = s in the 
tangent direction along the vortex line. 

Consider a small cubic box A(x) of size A > 6 and 
volume centred around the point x. We define the 
coarse-grained superfiuid vorticity field as 

a;,(x) = / s'd^, (9) 

This definition corresponds to the same coarse-graining 
procedure which was used in Ref [22] in a numerical cal- 
culation of a vortex tangle driven by an ABC normal 
fiow to show that the (coarse-grained) superfluid vorticity 
matches the vorticity of the normal flow. 

Note that ujs is nonzero only if the vortex lines are spa- 
tially organised. If the vortex lines point randomly in all 
directions, then each Cartesian component of is zero, 
because in each direction oriented vortex strands cancel 
each other out when summed algebraically. 

It is easy to check that the magnitude of is less than 
K times the local vortex line density L: 

I'^^WI = It^ / (10) 

^ iA(x) 

^ JA(x) ^ JA(x) 

because |s'| = 1. 

Prom the coarse-grained superfluid vorticity we can de- 
fine the (local) polarised vortex line density L||: 

AcL||(x) = |a;,(x)|, (11) 

Since 

L||(x)<L(x). (12) 

the missing part is a field which we call ix(x) and we 
have the (local) decomposition 

L(x) = Lx(x) + L||(x). (13) 

The smoothed field L|| filters the vortex tangle in k 
space, getting rid of short-wavelength Kelvin waves on the 
same vortex line. It also accounts for cancellation efiiects 
arising from vortex lines oriented in opposing directions. 
By construction, Ly is sufficiently organised that it defines 
the coarse-grained superfluid vorticity fleld u^g. Thus i|| 
reflects the superfluid velocity field in the inertial range 
at scales larger than A. At scales smaller than A the 
superfluid vorticity field Z/y is clearly not defined. 

Viccvcrsa, Lx = L — contributes to the vortex line 
density but not to the superfiuid vorticity and reflects 
the randomly oriented vortex lines. Note that we have 
not made any assumption about the relative amount of 



T^.0 
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wiggliness of Ly and ix- Notice that does not nec- 
essarily consist only of small loops (left over by vortex 
rcconncctions for example) or high energy Kelvin waves 
{k » A~^). Long filaments, provided they are randomly 
oriented with respect to their neighbours (so that they do 
not add up vortex length in the same direction), can be 
part of Lx- 

Passive vectors. The second feature of our model is 
the assumption that the unpolarised field Lx has some 
statistical properties of a classical passive vector field. It 
is well known [24] that passive vectors have a power spec- 
trum which obeys the /~^/^ law, and this is our expla- 
nation of the observed power spectrum of the vortex line 
density. 

At first it may seem contradictory to expect Lx to cor- 
respond to a active field while assuming a passive nature 
for We now argue that this active/passive distinction 
may result from a fundamental property of superfluid vor- 
tices which makes them different from classical vortices. 

In a classical fluid [25] the local time derivative of the 
vorticity is the combination of advection, (v • V)a;, and 
stretching, (u) • V)v. An example of the latter is the 
stretching of a classical vortex along its main axis which 
elongates it while squeezing it transversally, resulting in 
an increase of vorticity. In the superfluid vortex stretch- 
ing does not occur because the radius of the vortex core 
is fixed, determined by quantum mechanical constraints 
on the rotation. Superfluid vortices can become longer 
(for example, if energy is fed from the normal fluid, or, at 
T = 0, if the geometry changes keeping the total kinetic 
energy constant), but their core is rigid, thus they remain 
slender with respect to any typical turbulence scale. Su- 
perfluid vortex dynamics thus differs from the dynamics 
of classical vortices (for example see [25,26]). Neverthe- 
less, when superfluid vortices are assembled in polarised 
bundles, the classical vorticity enhancement that is de- 
scribed above can still be reproduced by stretching the 
whole bundle, as if vortices were material lines. For ex- 
ample, axisymmetric and non-axisymmetric oscillations of 
superfluid vortex bundles in the form of waves are known 
in the literature [27]. It is therefore reasonable to expect 
that the field w^, which result from polarised vortices, 
will mimic a classical active vorticity field at large enough 
scales. 

We now turn to the stretching of an unpolarised tan- 
gle by a large scale velocity field. If, as modelled above, 
superfluid vortices behave as material lines, it is also rea- 
sonable to expect that the tangle will remain unpolarised, 
and, due to the fluid incompressibility, the total length of 
lines will remain unchanged : in other words, Ly will be 
simply transported by the coarse-grained velocity. 

Viceversa, it is reasonable to assume that the high den- 
sity and large density fluctuations of L x have little impact 
on the dynamics of the polarised field Ly or Wg, that is to 
say that Lx does not advect Ly at scale large than A 
(where Ly is defined). This must be the case, because the 



velocity field induced by the unpolarised field must be very 
short-ranged, probably caused by multipolar sources 
with no contribution at first order 1/r. 

Consistency with meeisurements. — 

The spectrurn, of the polarised fi,eld. Our model is con- 
sistent only if we can show that the polarised vortex line 
density Ly gives a negligible contribution to the spectrum 
below IkHz, which is the observed second sound frequency 
range used by Roche et al. 

Following what said in the introduction, we approxi- 
mate the second sound spectrum P|| arising from as a 
white noise signal up to a viscous cutoff corresponding to 
few times the frequency of the Kolmogorov length scale r?: 



(14) 



The denominator is the full frequency span of an ideal 
second sound probe. Eq.l4 represents the (constant) 
power spectral density at all frequencies below the cut-off 
frequency. An ideal second sound probe is fixed in space. 
The smallest time scales which arc visible to the probe 
are produced by the smallest fiow structures (of typical 
size At]) which are advected past the probe at the local 
fluid velocity (approximated by the mean flow velocity). 
The inverse of this time scale gives the highest frequency 
of the signal seen by the ideal probe, which is indeed the 
frequency span. 

Let us find an upper bound for Py . In classical turbu- 
lence, both experimental [23] and numerical studies sug- 
gest that the vorticity u> roughly satisfies 

w2~2|IJ|^ (15) 

If we assume that this classical relation applies to the 
the normal fluid in the experiment [11], and if we make 
the further assumptions that normal fluid and superfluid 
are locked, » a;„, and that uLu = iVg, we have 



Thus 



Pi ~ ■=Lr, = 



V 



V 



(16) 



(17) 



where v — /i/p is the kinematic viscosity and /x the vis- 
cosity of helium. 

A bound for Py can be found by noticing that the total 
rate of dissipation of kinetic energy per unit volume in 
turbulent He II, pe, should be larger than the dissipation 
/zw^ which arises from the regular viscous dissipation in 
the normal fluid alone: 



where, since lOs ~ Wr? 



(18) 



(19) 
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Using the kinematic viscosity v = fi/ p, we find 



2^ 
V 



(20) 



In the next section, wc present different ways of evahiating 
the Koimogorov length -q &iT — l.QK and show that they 
all give the same order of magnitude for r/. For clarity, wc 
use the following classical expression (with the kinematic 
viscosity defined above): 



We have 



,1/4 



3 \ 1/4 



(21) 



(22) 



To evaluate this expression, we estimate the rate of 
kinetic energy dissipation e at the integral scale, « 
10~^m, for which Vq ~ Kms = 0.3m/s, and obtain 



T/3 



2.7m^ s" 



(23) 



At T = 1.6K, 

ll6 



8.9 X 10- 



s 1 and we obtain 



P|| < 4.3 X 10 °m~* s, which is much less then the observed 
spectral density P = 2 x lO^^m"^ s and is just above 
the instrumental noise level 0.5 x lOi^m"'' s (see Fig 4 of 
Ref. [11]). We conclude that the contribution of Ly to the 
observed spectrum is negligible. 

The Koimogorov length. We estimate the Koimogorov 
length from the expression ry = (i^*''/e)i/^. In principle wc 
can define three possible kinematic viscosities v* in our 
problem (all numerical values refer to T = 1.6K). The 
first is based on the total density p = p^-\- = 145kg m~^ 
and is V = pI p = S.'d X lO^^m^ s~i. The second is 
based on the normal fluid density p„ = 23.6kg m"*^ and 
is Vn = p/pn = 55 x 10~^m^ s^i. The third is the effi- 
cient kinematic viscosity j/' w 20x 10~^m^ s~i determined 
from towed grid experiments in turbulent He II, see fig. 8 
of Ref . [15] Using these three values, we obtain respec- 
tively T] « 0.7/Ltm, J] « 3/im and rj « 1pm. In all cases 
rj is of the order of magnitude of the intervortex spacing, 
S « Apm or at the most six times smaller. This is consis- 
tent with the argument of Vinen and Niemela [16] that 6 
and T] are likely to be of the same order of magnitude: the 
superfluid and the normal fluid are coupled throughout 
the inertial range. 

The energy. A key feature of our model is that the 

unpolariscd field L x gives a negligible contribution to the 
total energy of the flow compared to the polarised field . 
The estimate made in the introduction that the kinetic 
energy per unit volume is K ^ 6.5Jm^''' clearly refers to 
the polarised field Ly, because we obtained it using the 
condition u„ » ««: in this notation we rewrite 



Kn 



1 2 



1 2 



\{Pn 



6.5J m 



Ps J ' rms 
3 



(24) 



Let us estimate the energy contained in the unpolarised 
field ix • If we picture Lx as a random network of straight 
vortex lines, the energy per unit volume can be obtained 
following the integral procedure that leads to Eq.6 with 
b = 5: we multiply the length per unit volume times the 
integral of the square of the velocity field only up to a 
radial distance which is of the order of the intervortex 
spacing, because at this distance the total velocity fields 
of randomly oriented vortices cancel each other out. Using 
the numerical value provided in Eq.6 and Eq.2 for k. L, 
we get 

^2 



PsK 
477 



Lx In i6/a) « 0.06 Jm,-'' « K\ 



(25) 



The above integration procedure would not give the ki- 
netic energy per unit volume arising from the polarised 
field L||, because contributions of different vortex lines 
add up rather than cancel each other. This why we ex- 
pect K\\ » Kx to hold when a sufficient level of vortex 
polarisation is reached. 

The argument is made more clear if we consider a cylin- 
drical container of radius h and height h containing N 
straight vortex lines aligned along the axis. The number 
of vortex lines per unit area is L — N / (tt/i^) and the inter- 
vortex distance is5 = ^J-kH^ /N. Suppose that the vortices 
are oriented in the same direction (complete polarisation), 
forming a vortex bundle, as in a recent numerical calcula- 
tion [28] ; then the velocity fields of the vortices add up and 
create a total solid-body rotation velocity v = Or where 
il is obtained from 

Nk = i v-di= I VxvdS= / w-dS = 207^/l^ (26) 
J C J S J s 

which yields 



9. = NK/{2TTh^). 
The energy per unit volume is 



Ell 



16 



(27) 



(28) 



Now assume the opposite limit, that the vortices are 
randomly oriented in the positive or negative direction 
along the axis of the cylinder. The energy per unit volume 
is 



ln(l/(La2)), 



(We write it in terms of L rather than N because we want 
to take the limit of increasing the density rather than the 
number of vortices). Therefore 



Ex 



El 



^J^lnil/{La')) = ^^ln{5/a), 



(30) 



Clearly, the larger is the vortex line density L, the 
smaller is Ex with respect to E^^ (provided 5 > a of 
course). 
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Conclusion. — Our model of superfluid turbulence 
consists in dividing the tangle in a polarised field L\\ and 
a unpolarised field Lx such that the total vortex line den- 
sity is L = L|| + Lx. The polarised field makes up the 
(coarse-grained) superfluid vorticity field, such that the 
two fluids are coupled (w„ ~ in the inertial range ac- 
cording to current understanding of superfluid turbulence. 
The unpolarised field has some of the statistical properties 
of a passive vector field. What is observed with the sec- 
ond sound probe is mainly Lx, not L\\, which, as we have 
shown, brings a negligible contribution to the measured 
spectrum: 

Px»P\\. (31) 
However the polarised field Lu has more energy than L^: 



Ku » 



(32) 



That is why if we try to infer the vortex line density from 
Eq.8 we get an unrealistic high value. 

The model suggests the following picture of the turbu- 
lent tangle in the large scales (upper inertial range): large 
fiuctuations of the density of vortex line on top of a small 
modulation of polarisation which -nevertheless- controls 
the velocity field, energy cascade and passive advection of 
most of the vortex line density. 

We stress that the model which we propose is only an at- 
tempt to combine the information which is available from 
the experiment of Roche et al. [11] in a consistent sce- 
nario. If the model can be confirmed, the f~^^^ power law 
dependence of the fluctuations of the vortex line density 
reported in [11] should be considered as the inertial-range 
signature of the quantum nature of superfluid turbulence. 

Finally we remark that a decomposition of the vortex 
line density field in polarised and isotropic parts similar in 
spirit to what we have done, has been attempted by Lip- 
niacki [29]. Lipniacki's theory results in an Euler equation 
(motificd by the presence of friction) which is similar to 
Hall-Vinen equation for the macroscopic superfiuid veloc- 
ity [30] , which is coupled to a modified Vinen equation for 
the (more microscopic and isotropic) vortex line density 
L. 
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